3.112 \(\int \frac{1}{\sqrt{a+b x+c x^2} (d+e x+f x^2)} \, dx\)
Optimal. Leaf size=374 \[ \frac{\sqrt{2} f \tanh ^{-1}\left (\frac{4 a f+2 x \left (b f-c \left (\sqrt{e^2-4 d f}+e\right )\right )-b \left (\sqrt{e^2-4 d f}+e\right )}{2 \sqrt{2} \sqrt{a+b x+c x^2} \sqrt{2 a f^2+\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt{e^2-4 d f} \sqrt{2 a f^2+\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}-\frac{\sqrt{2} f \tanh ^{-1}\left (\frac{4 a f+2 x \left (b f-c \left (e-\sqrt{e^2-4 d f}\right )\right )-b \left (e-\sqrt{e^2-4 d f}\right )}{2 \sqrt{2} \sqrt{a+b x+c x^2} \sqrt{2 a f^2-\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt{e^2-4 d f} \sqrt{2 a f^2-\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}} \]
[Out]
-((Sqrt[2]*f*ArcTanh[(4*a*f - b*(e - Sqrt[e^2 - 4*d*f]) + 2*(b*f - c*(e - Sqrt[e^2 - 4*d*f]))*x)/(2*Sqrt[2]*Sq
rt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 - (c*e - b*f)*Sqrt[e^2 - 4*d*f]]*Sqrt[a + b*x + c*x^2])])/(Sqrt[e^2 - 4*d
*f]*Sqrt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 - (c*e - b*f)*Sqrt[e^2 - 4*d*f]])) + (Sqrt[2]*f*ArcTanh[(4*a*f - b*
(e + Sqrt[e^2 - 4*d*f]) + 2*(b*f - c*(e + Sqrt[e^2 - 4*d*f]))*x)/(2*Sqrt[2]*Sqrt[c*e^2 - 2*c*d*f - b*e*f + 2*a
*f^2 + (c*e - b*f)*Sqrt[e^2 - 4*d*f]]*Sqrt[a + b*x + c*x^2])])/(Sqrt[e^2 - 4*d*f]*Sqrt[c*e^2 - 2*c*d*f - b*e*f
+ 2*a*f^2 + (c*e - b*f)*Sqrt[e^2 - 4*d*f]])
________________________________________________________________________________________
Rubi [A] time = 0.57883, antiderivative size = 374, normalized size of antiderivative = 1.,
number of steps used = 5, number of rules used = 3, integrand size = 27, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.111, Rules used =
{983, 724, 206} \[ \frac{\sqrt{2} f \tanh ^{-1}\left (\frac{4 a f+2 x \left (b f-c \left (\sqrt{e^2-4 d f}+e\right )\right )-b \left (\sqrt{e^2-4 d f}+e\right )}{2 \sqrt{2} \sqrt{a+b x+c x^2} \sqrt{2 a f^2+\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt{e^2-4 d f} \sqrt{2 a f^2+\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}-\frac{\sqrt{2} f \tanh ^{-1}\left (\frac{4 a f+2 x \left (b f-c \left (e-\sqrt{e^2-4 d f}\right )\right )-b \left (e-\sqrt{e^2-4 d f}\right )}{2 \sqrt{2} \sqrt{a+b x+c x^2} \sqrt{2 a f^2-\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt{e^2-4 d f} \sqrt{2 a f^2-\sqrt{e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}} \]
Antiderivative was successfully verified.
[In]
Int[1/(Sqrt[a + b*x + c*x^2]*(d + e*x + f*x^2)),x]
[Out]
-((Sqrt[2]*f*ArcTanh[(4*a*f - b*(e - Sqrt[e^2 - 4*d*f]) + 2*(b*f - c*(e - Sqrt[e^2 - 4*d*f]))*x)/(2*Sqrt[2]*Sq
rt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 - (c*e - b*f)*Sqrt[e^2 - 4*d*f]]*Sqrt[a + b*x + c*x^2])])/(Sqrt[e^2 - 4*d
*f]*Sqrt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 - (c*e - b*f)*Sqrt[e^2 - 4*d*f]])) + (Sqrt[2]*f*ArcTanh[(4*a*f - b*
(e + Sqrt[e^2 - 4*d*f]) + 2*(b*f - c*(e + Sqrt[e^2 - 4*d*f]))*x)/(2*Sqrt[2]*Sqrt[c*e^2 - 2*c*d*f - b*e*f + 2*a
*f^2 + (c*e - b*f)*Sqrt[e^2 - 4*d*f]]*Sqrt[a + b*x + c*x^2])])/(Sqrt[e^2 - 4*d*f]*Sqrt[c*e^2 - 2*c*d*f - b*e*f
+ 2*a*f^2 + (c*e - b*f)*Sqrt[e^2 - 4*d*f]])
Rule 983
Int[1/(((a_) + (b_.)*(x_) + (c_.)*(x_)^2)*Sqrt[(d_.) + (e_.)*(x_) + (f_.)*(x_)^2]), x_Symbol] :> With[{q = Rt[
b^2 - 4*a*c, 2]}, Dist[(2*c)/q, Int[1/((b - q + 2*c*x)*Sqrt[d + e*x + f*x^2]), x], x] - Dist[(2*c)/q, Int[1/((
b + q + 2*c*x)*Sqrt[d + e*x + f*x^2]), x], x]] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[e
^2 - 4*d*f, 0] && NeQ[c*e - b*f, 0] && PosQ[b^2 - 4*a*c]
Rule 724
Int[1/(((d_.) + (e_.)*(x_))*Sqrt[(a_.) + (b_.)*(x_) + (c_.)*(x_)^2]), x_Symbol] :> Dist[-2, Subst[Int[1/(4*c*d
^2 - 4*b*d*e + 4*a*e^2 - x^2), x], x, (2*a*e - b*d - (2*c*d - b*e)*x)/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a,
b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[2*c*d - b*e, 0]
Rule 206
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1*ArcTanh[(Rt[-b, 2]*x)/Rt[a, 2]])/(Rt[a, 2]*Rt[-b, 2]), x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])
Rubi steps
\begin{align*} \int \frac{1}{\sqrt{a+b x+c x^2} \left (d+e x+f x^2\right )} \, dx &=\frac{(2 f) \int \frac{1}{\left (e-\sqrt{e^2-4 d f}+2 f x\right ) \sqrt{a+b x+c x^2}} \, dx}{\sqrt{e^2-4 d f}}-\frac{(2 f) \int \frac{1}{\left (e+\sqrt{e^2-4 d f}+2 f x\right ) \sqrt{a+b x+c x^2}} \, dx}{\sqrt{e^2-4 d f}}\\ &=-\frac{(4 f) \operatorname{Subst}\left (\int \frac{1}{16 a f^2-8 b f \left (e-\sqrt{e^2-4 d f}\right )+4 c \left (e-\sqrt{e^2-4 d f}\right )^2-x^2} \, dx,x,\frac{4 a f-b \left (e-\sqrt{e^2-4 d f}\right )-\left (-2 b f+2 c \left (e-\sqrt{e^2-4 d f}\right )\right ) x}{\sqrt{a+b x+c x^2}}\right )}{\sqrt{e^2-4 d f}}+\frac{(4 f) \operatorname{Subst}\left (\int \frac{1}{16 a f^2-8 b f \left (e+\sqrt{e^2-4 d f}\right )+4 c \left (e+\sqrt{e^2-4 d f}\right )^2-x^2} \, dx,x,\frac{4 a f-b \left (e+\sqrt{e^2-4 d f}\right )-\left (-2 b f+2 c \left (e+\sqrt{e^2-4 d f}\right )\right ) x}{\sqrt{a+b x+c x^2}}\right )}{\sqrt{e^2-4 d f}}\\ &=-\frac{\sqrt{2} f \tanh ^{-1}\left (\frac{4 a f-b \left (e-\sqrt{e^2-4 d f}\right )+2 \left (b f-c \left (e-\sqrt{e^2-4 d f}\right )\right ) x}{2 \sqrt{2} \sqrt{c e^2-2 c d f-b e f+2 a f^2-(c e-b f) \sqrt{e^2-4 d f}} \sqrt{a+b x+c x^2}}\right )}{\sqrt{e^2-4 d f} \sqrt{c e^2-2 c d f-b e f+2 a f^2-(c e-b f) \sqrt{e^2-4 d f}}}+\frac{\sqrt{2} f \tanh ^{-1}\left (\frac{4 a f-b \left (e+\sqrt{e^2-4 d f}\right )+2 \left (b f-c \left (e+\sqrt{e^2-4 d f}\right )\right ) x}{2 \sqrt{2} \sqrt{c e^2-2 c d f-b e f+2 a f^2+(c e-b f) \sqrt{e^2-4 d f}} \sqrt{a+b x+c x^2}}\right )}{\sqrt{e^2-4 d f} \sqrt{c e^2-2 c d f-b e f+2 a f^2+(c e-b f) \sqrt{e^2-4 d f}}}\\ \end{align*}
Mathematica [A] time = 1.50838, size = 376, normalized size = 1.01 \[ \frac{\sqrt{2} f \left (\frac{\tanh ^{-1}\left (\frac{4 a f-b \left (\sqrt{e^2-4 d f}+e-2 f x\right )-2 c x \left (\sqrt{e^2-4 d f}+e\right )}{2 \sqrt{2} \sqrt{a+x (b+c x)} \sqrt{f \left (2 a f-b \left (\sqrt{e^2-4 d f}+e\right )\right )+c \left (e \sqrt{e^2-4 d f}-2 d f+e^2\right )}}\right )}{\sqrt{f \left (2 a f-b \left (\sqrt{e^2-4 d f}+e\right )\right )+c \left (e \sqrt{e^2-4 d f}-2 d f+e^2\right )}}-\frac{\tanh ^{-1}\left (\frac{4 a f+b \left (\sqrt{e^2-4 d f}-e+2 f x\right )+2 c x \left (\sqrt{e^2-4 d f}-e\right )}{2 \sqrt{2} \sqrt{a+x (b+c x)} \sqrt{f \left (2 a f+b \sqrt{e^2-4 d f}+b (-e)\right )+c \left (-e \sqrt{e^2-4 d f}-2 d f+e^2\right )}}\right )}{\sqrt{f \left (2 a f+b \sqrt{e^2-4 d f}+b (-e)\right )+c \left (-e \sqrt{e^2-4 d f}-2 d f+e^2\right )}}\right )}{\sqrt{e^2-4 d f}} \]
Antiderivative was successfully verified.
[In]
Integrate[1/(Sqrt[a + b*x + c*x^2]*(d + e*x + f*x^2)),x]
[Out]
(Sqrt[2]*f*(ArcTanh[(4*a*f - 2*c*(e + Sqrt[e^2 - 4*d*f])*x - b*(e + Sqrt[e^2 - 4*d*f] - 2*f*x))/(2*Sqrt[2]*Sqr
t[c*(e^2 - 2*d*f + e*Sqrt[e^2 - 4*d*f]) + f*(2*a*f - b*(e + Sqrt[e^2 - 4*d*f]))]*Sqrt[a + x*(b + c*x)])]/Sqrt[
c*(e^2 - 2*d*f + e*Sqrt[e^2 - 4*d*f]) + f*(2*a*f - b*(e + Sqrt[e^2 - 4*d*f]))] - ArcTanh[(4*a*f + 2*c*(-e + Sq
rt[e^2 - 4*d*f])*x + b*(-e + Sqrt[e^2 - 4*d*f] + 2*f*x))/(2*Sqrt[2]*Sqrt[f*(-(b*e) + 2*a*f + b*Sqrt[e^2 - 4*d*
f]) + c*(e^2 - 2*d*f - e*Sqrt[e^2 - 4*d*f])]*Sqrt[a + x*(b + c*x)])]/Sqrt[f*(-(b*e) + 2*a*f + b*Sqrt[e^2 - 4*d
*f]) + c*(e^2 - 2*d*f - e*Sqrt[e^2 - 4*d*f])]))/Sqrt[e^2 - 4*d*f]
________________________________________________________________________________________
Maple [B] time = 0.359, size = 761, normalized size = 2. \begin{align*} \text{result too large to display} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(1/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x)
[Out]
-1/(-4*d*f+e^2)^(1/2)*2^(1/2)/(((-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-b*e*f-2*c*d*f+c*e^2)/f^2
)^(1/2)*ln((((-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-b*e*f-2*c*d*f+c*e^2)/f^2+(c*(-4*d*f+e^2)^(1
/2)+b*f-c*e)/f*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)+1/2*2^(1/2)*(((-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e+2
*a*f^2-b*e*f-2*c*d*f+c*e^2)/f^2)^(1/2)*(4*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)^2*c+4*(c*(-4*d*f+e^2)^(1/2)+b*f-c*
e)/f*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)+2*((-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-b*e*f-2*c*d*f+
c*e^2)/f^2)^(1/2))/(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f))+1/(-4*d*f+e^2)^(1/2)*2^(1/2)/((-(-4*d*f+e^2)^(1/2)*b*f+(
-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-b*e*f-2*c*d*f+c*e^2)/f^2)^(1/2)*ln(((-(-4*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*
c*e+2*a*f^2-b*e*f-2*c*d*f+c*e^2)/f^2+1/f*(-c*(-4*d*f+e^2)^(1/2)+b*f-c*e)*(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f)+1/2*
2^(1/2)*((-(-4*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-b*e*f-2*c*d*f+c*e^2)/f^2)^(1/2)*(4*(x+1/2*(e+
(-4*d*f+e^2)^(1/2))/f)^2*c+4/f*(-c*(-4*d*f+e^2)^(1/2)+b*f-c*e)*(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f)+2*(-(-4*d*f+e^
2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-b*e*f-2*c*d*f+c*e^2)/f^2)^(1/2))/(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f))
________________________________________________________________________________________
Maxima [F(-2)] time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="maxima")
[Out]
Exception raised: ValueError
________________________________________________________________________________________
Fricas [B] time = 45.0179, size = 22873, normalized size = 61.16 \begin{align*} \text{result too large to display} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="fricas")
[Out]
1/4*sqrt(2)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f + (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*
b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt
((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5
+ (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^
4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*
a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3
- (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (
b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*
f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e
^2)*f))*log((2*(b^2*d - a*b*e)*f^2 + sqrt(2)*(c^2*d*e^3 - 4*a*b*d*f^3 + (4*b*c*d^2 + 4*a*c*d*e + a*b*e^2)*f^2
- (4*c^2*d^2*e + b*c*d*e^2 + a*c*e^3)*f - (c^3*d^3*e^3 - b*c^2*d^2*e^4 + a*c^2*d*e^5 + 4*(2*a^2*b*d^2 - a^3*d*
e)*f^4 + (2*a^2*b*d*e^2 + a^3*e^3 + 8*(b^3 - 2*a*b*c)*d^3 - 4*(3*a*b^2 - a^2*c)*d^2*e)*f^3 + (8*b*c^2*d^4 - a^
2*b*e^4 - 4*(3*b^2*c - a*c^2)*d^3*e - 2*(b^3 - 10*a*b*c)*d^2*e^2 + (3*a*b^2 - 5*a^2*c)*d*e^3)*f^2 - (4*c^3*d^4
*e - 2*b*c^2*d^3*e^2 + 4*a*b*c*d*e^4 - a^2*c*e^5 - (3*b^2*c - 5*a*c^2)*d^2*e^3)*f)*sqrt((c^2*e^2 - 2*b*c*e*f +
b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d
^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c
^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c -
a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*
e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*
e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))*sqrt(c*x^2 + b*x + a)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f + (c^2
*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 -
4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^
3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8
*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e
+ (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c
+ 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d
^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4
)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4
*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)) - 2*(b*c*d*e - a*c*e^2)*f + ((4*b*c*d - b^2*e)*f^2
- (4*c^2*d*e - b*c*e^2)*f)*x - (8*a^3*d*f^4 - 2*(4*a^2*b*d*e + a^3*e^2 - 4*(a*b^2 - 2*a^2*c)*d^2)*f^3 + 2*(4*
a*c^2*d^3 - 4*a*b*c*d^2*e + a^2*b*e^3 - (a*b^2 - 6*a^2*c)*d*e^2)*f^2 - 2*(a*c^2*d^2*e^2 - a*b*c*d*e^3 + a^2*c*
e^4)*f + (4*a^2*b*d*f^4 - (4*a*b^2*d*e + a^2*b*e^2 - 4*(b^3 - 2*a*b*c)*d^2)*f^3 + (4*b*c^2*d^3 - 4*b^2*c*d^2*e
+ a*b^2*e^3 - (b^3 - 6*a*b*c)*d*e^2)*f^2 - (b*c^2*d^2*e^2 - b^2*c*d*e^3 + a*b*c*e^4)*f)*x)*sqrt((c^2*e^2 - 2*
b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2
*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c
+ 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8
*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 +
2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*
c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/x) - 1/4*sqrt(2)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f +
(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^
3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^
3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2
- 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d
^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b
^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c
^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d
*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2
- (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f))*log((2*(b^2*d - a*b*e)*f^2 - sqrt(2)*(c^2*d*e^
3 - 4*a*b*d*f^3 + (4*b*c*d^2 + 4*a*c*d*e + a*b*e^2)*f^2 - (4*c^2*d^2*e + b*c*d*e^2 + a*c*e^3)*f - (c^3*d^3*e^3
- b*c^2*d^2*e^4 + a*c^2*d*e^5 + 4*(2*a^2*b*d^2 - a^3*d*e)*f^4 + (2*a^2*b*d*e^2 + a^3*e^3 + 8*(b^3 - 2*a*b*c)*
d^3 - 4*(3*a*b^2 - a^2*c)*d^2*e)*f^3 + (8*b*c^2*d^4 - a^2*b*e^4 - 4*(3*b^2*c - a*c^2)*d^3*e - 2*(b^3 - 10*a*b*
c)*d^2*e^2 + (3*a*b^2 - 5*a^2*c)*d*e^3)*f^2 - (4*c^3*d^4*e - 2*b*c^2*d^3*e^2 + 4*a*b*c*d*e^4 - a^2*c*e^5 - (3*
b^2*c - 5*a*c^2)*d^2*e^3)*f)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d
*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c
)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3
*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*
e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5
+ (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))*sqrt(c*x^2 +
b*x + a)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f + (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*
d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((
c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 +
(b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4
- 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*
c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 -
(a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^
3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^
3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2
)*f)) - 2*(b*c*d*e - a*c*e^2)*f + ((4*b*c*d - b^2*e)*f^2 - (4*c^2*d*e - b*c*e^2)*f)*x - (8*a^3*d*f^4 - 2*(4*a^
2*b*d*e + a^3*e^2 - 4*(a*b^2 - 2*a^2*c)*d^2)*f^3 + 2*(4*a*c^2*d^3 - 4*a*b*c*d^2*e + a^2*b*e^3 - (a*b^2 - 6*a^2
*c)*d*e^2)*f^2 - 2*(a*c^2*d^2*e^2 - a*b*c*d*e^3 + a^2*c*e^4)*f + (4*a^2*b*d*f^4 - (4*a*b^2*d*e + a^2*b*e^2 - 4
*(b^3 - 2*a*b*c)*d^2)*f^3 + (4*b*c^2*d^3 - 4*b^2*c*d^2*e + a*b^2*e^3 - (b^3 - 6*a*b*c)*d*e^2)*f^2 - (b*c^2*d^2
*e^2 - b^2*c*d*e^3 + a*b*c*e^4)*f)*x)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*
a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2
- 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^
2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*
c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2
*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/x)
+ 1/4*sqrt(2)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f - (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*
a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sq
rt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^
5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(
b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 -
2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e
^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 +
(b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*
d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d
*e^2)*f))*log((2*(b^2*d - a*b*e)*f^2 + sqrt(2)*(c^2*d*e^3 - 4*a*b*d*f^3 + (4*b*c*d^2 + 4*a*c*d*e + a*b*e^2)*f^
2 - (4*c^2*d^2*e + b*c*d*e^2 + a*c*e^3)*f + (c^3*d^3*e^3 - b*c^2*d^2*e^4 + a*c^2*d*e^5 + 4*(2*a^2*b*d^2 - a^3*
d*e)*f^4 + (2*a^2*b*d*e^2 + a^3*e^3 + 8*(b^3 - 2*a*b*c)*d^3 - 4*(3*a*b^2 - a^2*c)*d^2*e)*f^3 + (8*b*c^2*d^4 -
a^2*b*e^4 - 4*(3*b^2*c - a*c^2)*d^3*e - 2*(b^3 - 10*a*b*c)*d^2*e^2 + (3*a*b^2 - 5*a^2*c)*d*e^3)*f^2 - (4*c^3*d
^4*e - 2*b*c^2*d^3*e^2 + 4*a*b*c*d*e^4 - a^2*c*e^5 - (3*b^2*c - 5*a*c^2)*d^2*e^3)*f)*sqrt((c^2*e^2 - 2*b*c*e*f
+ b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)
*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2
*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c
- a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c
)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^
2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))*sqrt(c*x^2 + b*x + a)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f - (c
^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3
- 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*
d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 -
8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2
*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2
*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3
*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e
^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 -
(4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)) - 2*(b*c*d*e - a*c*e^2)*f + ((4*b*c*d - b^2*e)*f
^2 - (4*c^2*d*e - b*c*e^2)*f)*x + (8*a^3*d*f^4 - 2*(4*a^2*b*d*e + a^3*e^2 - 4*(a*b^2 - 2*a^2*c)*d^2)*f^3 + 2*(
4*a*c^2*d^3 - 4*a*b*c*d^2*e + a^2*b*e^3 - (a*b^2 - 6*a^2*c)*d*e^2)*f^2 - 2*(a*c^2*d^2*e^2 - a*b*c*d*e^3 + a^2*
c*e^4)*f + (4*a^2*b*d*f^4 - (4*a*b^2*d*e + a^2*b*e^2 - 4*(b^3 - 2*a*b*c)*d^2)*f^3 + (4*b*c^2*d^3 - 4*b^2*c*d^2
*e + a*b^2*e^3 - (b^3 - 6*a*b*c)*d*e^2)*f^2 - (b*c^2*d^2*e^2 - b^2*c*d*e^3 + a*b*c*e^4)*f)*x)*sqrt((c^2*e^2 -
2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 +
2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*
c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 -
8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2
+ 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*
b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/x) - 1/4*sqrt(2)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f
- (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*
d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*
c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e
^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)
*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a
*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b
*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)
*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^
2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f))*log((2*(b^2*d - a*b*e)*f^2 - sqrt(2)*(c^2*d*
e^3 - 4*a*b*d*f^3 + (4*b*c*d^2 + 4*a*c*d*e + a*b*e^2)*f^2 - (4*c^2*d^2*e + b*c*d*e^2 + a*c*e^3)*f + (c^3*d^3*e
^3 - b*c^2*d^2*e^4 + a*c^2*d*e^5 + 4*(2*a^2*b*d^2 - a^3*d*e)*f^4 + (2*a^2*b*d*e^2 + a^3*e^3 + 8*(b^3 - 2*a*b*c
)*d^3 - 4*(3*a*b^2 - a^2*c)*d^2*e)*f^3 + (8*b*c^2*d^4 - a^2*b*e^4 - 4*(3*b^2*c - a*c^2)*d^3*e - 2*(b^3 - 10*a*
b*c)*d^2*e^2 + (3*a*b^2 - 5*a^2*c)*d*e^3)*f^2 - (4*c^3*d^4*e - 2*b*c^2*d^3*e^2 + 4*a*b*c*d*e^4 - a^2*c*e^5 - (
3*b^2*c - 5*a*c^2)*d^2*e^3)*f)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2
*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3
*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a
^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^
2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^
5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))*sqrt(c*x^2
+ b*x + a)*sqrt((c*e^2 + 2*a*f^2 - (2*c*d + b*e)*f - (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*
b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt
((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5
+ (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^
4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*
a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3
- (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (
b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*
f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e
^2)*f)) - 2*(b*c*d*e - a*c*e^2)*f + ((4*b*c*d - b^2*e)*f^2 - (4*c^2*d*e - b*c*e^2)*f)*x + (8*a^3*d*f^4 - 2*(4*
a^2*b*d*e + a^3*e^2 - 4*(a*b^2 - 2*a^2*c)*d^2)*f^3 + 2*(4*a*c^2*d^3 - 4*a*b*c*d^2*e + a^2*b*e^3 - (a*b^2 - 6*a
^2*c)*d*e^2)*f^2 - 2*(a*c^2*d^2*e^2 - a*b*c*d*e^3 + a^2*c*e^4)*f + (4*a^2*b*d*f^4 - (4*a*b^2*d*e + a^2*b*e^2 -
4*(b^3 - 2*a*b*c)*d^2)*f^3 + (4*b*c^2*d^3 - 4*b^2*c*d^2*e + a*b^2*e^3 - (b^3 - 6*a*b*c)*d*e^2)*f^2 - (b*c^2*d
^2*e^2 - b^2*c*d*e^3 + a*b*c*e^4)*f)*x)*sqrt((c^2*e^2 - 2*b*c*e*f + b^2*f^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 -
2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^
2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*
b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^
2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a
^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/x
)
________________________________________________________________________________________
Sympy [F] time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{1}{\sqrt{a + b x + c x^{2}} \left (d + e x + f x^{2}\right )}\, dx \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c*x**2+b*x+a)**(1/2)/(f*x**2+e*x+d),x)
[Out]
Integral(1/(sqrt(a + b*x + c*x**2)*(d + e*x + f*x**2)), x)
________________________________________________________________________________________
Giac [F(-1)] time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="giac")
[Out]
Timed out